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Abstract 

Within the Lorentz invariant formalizm for description of neutrino evolution 
in electromagnetic fields and matter we consider neutrino spin oscillations in the 
electromagnetic wave with varying amplitude and in "castle wall" magnetic field. 
It is shown for the first time that the parametric resonances of neutrino oscillations 
in such systems can occur. 
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Experimental studies of solar, atmospheric and reactor neutrinos over the past several 
years provide almost certain indications that neutrinos oscillate, have masses and mix. 
These new properties of neutrinos, if confirmed by better statistics in the proceeding and 
forthcomming experiments, will require a significant departure from the Standard Model. 
It is well known [1] that massive neutrinos can have nonvanishing magnetic moment. For 
example, in the Standard Model supplied with SU(2)-smg\et right-handed neutrino the 
one-loop radiative correction generates neutrino magnetic moment which is proportional 
to neutrino mass. There are plenty of models [2] which predict much large magnetic 
moment for neutrinos. 

In the light of new developments in neutrino physics, to get better understanding of 
the electromagnetic properties of neutrinos is an important task. In general, there are two 
aspects of this problem. The first one is connected with values of neutrino electromagnetic 
form - factors. The second aspect of the problem implies consideration of influence of 
external electromagnetic fields, which could be presented in various environments, on 
neutrino possessing non - vanishing electromagnetic moments. The most important of 
the latter are the magnetic and electric dipole moments. 

The majority of the previously performed studies of neutrino conversions and oscilla- 
tions in electromagnetic fields deal with the case of constant transversal magnetic field or 
constant transversal twisting magnetic fields (see references [1-18] of paper [3]). Recently 
we have developed [3-5] the Loretz invariant formalism for description of neutrino spin 
evolution that enables one to consider neutrino oscillations in the presence of an arbi- 
trary electromagnetic fields. Within the proposed approach it becomes possible to study 
neutrino spin evolution in an electromagnetic wave and the new types of resonance in 
the neutrino oscillations in the wave field and in some other combinations of fields have 
been predicted. The latest application [6] of this formalism enables us to consider neu- 
trino oscillations in electromagnetic fields in moving and polarised matter and show that 
matter effects on neutrino oscillations in the case of relativistic motion of matter could 
sufficiently depend on matter motion. 

The aim of this paper is to continue the study of the neutrino spin oscillations in 
different electromagnetic fields and consider the cases of electromagnetic wave with peri- 
odically varying amplitude (see also [7]) and periodic "castle - wall" magnetic field. It will 
be shown that under certain conditions the enhancement of neutrino oscillations could be 
achieved due to the parametric resonance effect. 

The possibility of the parametric resonance of neutrino flavour oscillations in mat- 
ter with periodic variation of density was considered previously in [8]. It should be also 



pointed here that the conditions for a total neutrino flavour conversion in a medium 
consisting of two or three constant density layers were derived in [9] . For the recent dis- 
cussion on the physical interpretations of these two mechanisms of increasing of neutrino 
conversion see refs.[9-ll]. 

First we consider the case when the case of electromagnetic wave, with the amplitude 
of which being not constant but periodically varying in time, and show (see also ref. [7]) 
that under certain conditions the parametric resonance of neutrino oscillations can occur 
in such a system. Note that in order to investigate this phenomenon we have to use the 
Lorentz invariant formalisms for neutrino evolution that have been developed in [3-5]. 

Let us consider evolution of a system v = (is + ,is_) composed of two neutrinos of 
different helicities in presence of a field of circular polarised electromagnetic wave with 
periodically varying amplitude. Here neutrinos u + and i/_ correspond to positive and 
negative helicity states that are determined by projection operators 

In general case it is important to distinguish helicity states v = (z/ + ,z/_) and chirality 
states v = (u L , u R ). The latter are determined by the projection operators Pl,r = (1 =F 
75)/2. The evolution of such a system is given [3-5] by the following Schrodinger type 
equation 

dv 

i— = Hu, H = pa 3 + E(t)(ai cos?/' — o"2 sin tp), (2) 

p = -^ + ^^, E{t) =»B{t){l- (3 cos 4>). 

Here the three parameters, A = A(9) being a function of vacuum mixing angle, V = 
V(n e ff) being the difference of neutrino effective potentials in matter, and Am 2 being the 
neutrino masses squired difference depend on the considered type of neutrino conversion 
process. The electromagnetic field is determined in the laboratory frame of reference by its 
frequency, u, the phase at the point of the neutrino location, ip = gut(l — (3/ (3q cos <j)),(g = 
±1), and the amplitude, B(t), which is a function of time. The wave speed in matter 
could be less than the speed of light in vacuum (f3 < 1), and is the angle between 
the neutrino speed /3 and the direction of the wave propagation. In the derivation of 
the Hamiltonian of eq.(2) terms proportional to -^ <S 1, 7 = (1 — /J 2 ) -1 / 2 , and also an 
oscillating function of time in the diagonal part are neglected (for details see [3-5]). 



In order to study phenomenon of the parametric resonance of neutrino spin oscillations 
in such a wave we suppose that the amplitude B(t) is given by 

B(t) = B(l + hf(t)), (|/i|«l), (3) 

where /(£) is an arbitrary function of time and h is a small dimensionless quantity of 
not fixed sign. It is convenient to introduce the evolution operator which determines the 
neutrino state at time t 

u(t) = U(t)u(0) 

if the initial neutrino state was v(0). Using the Hamiltonian of eq.(2) we get the following 
equation for the evolution operator: 

t/(t) = t[-pa 3 + (-E + e/(t))((Ti cos^ - a 2 sm^)]U(t), (4) 

where ce = —EqH. In analogy with the case of the electromagnetic wave with non - 
varying amplitude [3,4] the solution of eq.(4) can be written in the form 

U(t)=U, 3 (t)U l {t)F(t), (5) 

where Ug 3 (t) = exp(ia 3 ^) is the rotation operator around the axis e 3 which is parallel 
with the direction of the neutrino propagation, and UAt) = exp(ialt) is the rotation 
operator around the vector I = (—E ,0,p — ^). It should be noted that the solution of 
eq.(4) for the case of constant amplitude of the wave field (e = 0) is given by the operator 
U (t)=U gs (t)U J (t). 

From (4) and (5) it follows that the equation for the operator F(t) is 

F(t)=ieH e (t)F(t), (6) 

where 

H £ (t) = (ay(t))f(t), 

y= (l-2\ 2 3 sm 2 ttt, A 3 sin2fit, 2AiA 3 sin 2 fit), 

and the unit vector A is given by its components in the unit orthogonal basis A = Aiei + 
A 2 e 2 + A 3 e 3 = ^, Q = \l\. The detailed analysis of evolution of the solution of eq.(6) can 
be found in [12]. Here we comment only on the main steps. Using the smallness of e we 
expand the solution of eq.(6) in powers of this parameter: 

oo 

F = Y / e k F^ k \ (7) 

fc=0 



where F^- ' — 1 is a unit matrix. For operators F^ k ' the recurrent formula is straightfor- 
ward: 

t 

F (fe+1) (t) = i f H £ (T)F (k) (T)dT. (8) 



Skipping further technical details to the first order in e we get 

F(t) = 1 + ie(ax(t)) + 0(e 2 ), (9) 

where 

t 

x(t)=Jy(r)f(r)dr. 
o 
Thus for the probability of neutrino conversion z/j_ <->■ i>j + we get 

Pij = | < v+\Ug 3 (t)Ui{t)F(t)\v- > | 2 = (10) 

A 2 sin 2 fit + 2eAi (xi (i) cos fit + A3X2 (t) sin Qt) sin fit. 
Note that fi is nothing but the mean osillation frequency of the neutrino system. 

For further evaluation of solution of eq.(6) we have to specify the form of the function 
/(t). As it is mentioned above, the purpose of this study is to examine the case when 
the parametric resonance for neutrino oscillations in electromagnetic wave with varying 
amplitude could appear. Having in mind a simple analogy ( see, for example, [10]) be- 
tween oscillations in neutrino system and oscillations of a classical pendulum [13], it is 
reasonable to suppose that the principal parametric resonance appears when the ampli- 
tude modulation function f(t) is oscillating in time with frequency approximately equal 
to the twice mean oscillation frequency of the system. That is why we choose the function 
fit) to be 

/(t) = sin2fit (11) 

and get for the neutrino oscillations probability 

f\ A 2 

Pij = [X 2 1 + s\ 1 \lt+— i(l-^)sin2fit]sin 2 fit. (12) 

It follows that in the case Ai£ > the second term increases with increase of time t, so 
that the amplitude of the neutrino conversion probability may becomes close to unity. 
This is the effect of the parametric resonance in neutrino system in the electromagnetic 
wave with modulated amplitude that may enhance neutrino oscillation amplitude even for 
rather small values of the neutrino magnetic moment \x and strength of the electromagnetic 



field and also when the system is far away from the region of ordinary spin (or spin-flavour) 
neutrino resonance. 

Let us consider the case when parameters of the neutrino system are far beyond the 
region of the ordinary neutrino resonance. Supposing that \li\ = 0. 1 1 Z3 1 . In this case 
the maximal neutrino conversion probability, for non - varying (h = 0) amplitude of the 
electromagnetic field, is indeed small, 

72 

P ljmax (h = 0) = w -^ < 1. (13) 

n + '3 

If the amplitude of the field is modulated in accordance with eq.(ll) the estimation for the 

critical time t cr for which the probability P^ could become close to unity for an arbitrary 

values of the neutrino mixing angle, masses, energy, and density of matter gives 

to- « • (14) 

en 1 
It means that the parametric resonance enhancement of the neutrino oscillations occurs 
after neutrinos travel a distance under the influence of electromagnetic wave with mod- 
ulated amplitude. For the further numerical estimates me chose the strength of the field 
B rs 10 3 G, the value of the neutrino magnetic moment fi ~ 10~ 10 /i,B (/Ab is the Borh 
magneton), and also suppose that neutrino is propagating in the direction opposite to one 
of the electromagnetic wave. Thus, from (14) we obtain, for the small time variations of 
the field amplitude (taking, for instance \h\ = 0.1), that the conversion probability could 
attain the value Pij ~ 1 after the neutrino has travelled a distance 

Xrnin = 10~ U m. (15) 

Let us consider now a system of two helicity neutrino states with mixing v = (z/ + , i/_). 
In general the neutrino components can belong to different flavours. The evolution of 
such a system in matter under the influence of transversal magnetic field is described by 
the Schrodinger equation: 

«! = J7„ (16) 



with Hamiltonian 



H = -nB{t)a, + pa 3 , p = ^-A - V -f- (17) 



Now we are going to discuss a possibility of parametric resonance amplification that could 
exist under the influence of constant in time but periodically varying in space magnetic 
field. We consider a special case when the strength of the field varies along the neutrino 
path and the dependence on the neutrino coordinate t is given by, 

1 B 2 , 7\ < t < Tx + T 2 y J 

here B\ and B 2 are constants. The field strength is periodic "castle wall" function with 
period T = Ti + T 2 : B(t + T) = B(t). The effective Hamiltonian is also a periodic function 
of time with the same period T: 

H(t + T) = H(t), (19) 

and 

Hi, H 2 are the constant operators. It should be noted that the similar approach to 
neutrino oscillations in the case of periodic step - function ("castle wall") density was 
described in details in the first paper of ref. [10] 

Let us define the evolution operators for the two time intervals (0, Ti) and (Ti,T): 

U a = exp(-iH a T a ), a = 1,2. (21) 

Then the operator giving the neutrino evolution for the whole period T is 

U T = U 2 Ui. (22) 

For the further consideration it is convenient to introduce the unit vectors 

n a = (£ o ,O,-p) = (sin20 a ,O,-cos20 o ), a = 1,2, (23) 

UJn. 



where E a = \xB a , u a = Jp 2 + E 2 and 6 a are the effective neutrino mixing angles account- 
ing for interaction with matter and magnetic field. Then for effective Hamiltonians we 
can write 

H a = u a {a-n a ), (24) 



As it is easy to see, the eigenvalues of H a are ±u a . Using eqs. (18-24) we get the 
followig expressions for the neutrino evolution operator over the whole period T: 



where 



and 



U T = Y-i(a-X) = exp(-i(a-n x )$), (25) 

Y = cic 2 - (Hi ■ n 2 )sis 2 , (26) 

— * 

X = sic 2 ni + s 2 c x H 2 - (Hi x n 2 )sis 2 , (27) 

$ = arcsinF = arccosX, (28) 

n x = |, X=\X\. (29) 

s a = sin0 a , c a = cos0 a , (p a = uj a T a , a =1,2. (30) 

From (26), (27) it follows that 

HiH 2 = —(BiB 2 + p 2 ) =cos2(9i-0 2 ), (31) 

LOiUJ 2 

Hi xH 2 = ^(0,p(Bi-B 2 ),0) = (0,sm2($i-e 2 ),0). (32) 

UJlUJ 2 

It should be noted also that Y 2 + X 2 = 1 due to unitarity of the operator Ut- The vector 
X can be expressed in its components as 

X = (- S -^^ - S ^^,~P—(E 2 - Ei),p(— + ^)). (33) 

V OJl <jJ 2 UJ\UJ 2 UJ\ UJ 2 ' 

Let us consider the case when neutrino path in matter and periodic "castle wall" 
magnetic field is exactly equal to a whole number of periods, so that t = nT, n e N. The 
neutrino evolution operator for n periods is equal to Ut in the nth power. Assume that 
initially the two - level neutrino system concerned occupies the state z/_. 

U nT = exp(-i(a ■ H x )n®). (34) 

Assume that initially (t = 0) the two - level neutrino system involved occupies the 
state V-. The transition probability to the state v + after a time t > is determined by 
the evolution operator Ut' 



P(t) = I <iy + \U(t)\u_ > | 2 . (35) 

For the case t = nT using eqs.(34),(35) we get 



P (nT) = ^r^T^l ^ 2 (n$) = v^W^V2 ^ 2 ($-). (36) 



Xl + X 2 . 2 X-, + X 2 . 2 £ 

x i + xi + xi * n (n4>) = xFTxFTxi sm ( *r 

This expression coincides with one obtained in the first paper of ref.[10] for neutrino 
oscillation probability in step - function profile of matter density. However, it is obvious 
that in our case Xj and $ have different meaning. 

Expression (36) also looks like the neutrino oscillations probability for the case of not 
varying in space (non-periodic) magnetic field. However, there is an important distinction 
from the case of non-periodic magnetic field. The pre-sine factor in the latter case equals 
g 2+ -2 which is small if the value fxB is rather small. Contrary to the case of non-periodic 
magnetic field, in the case of B\ ^ B 2 the pre-sine factor need not be small even if sin 2 29 a 
are small. For particular choice of parameters the pre-sine term in eq.(36) becomes equal 
to unity. It is this case which we call the parametric enhancement of neutrino oscillations 
in periodic "castle wall" magnetic field. The resonance condition is 

Xl = p 2 (^ + ^f = 0, (37) 

Now let us consider the case when the neutrino oscillations can be considerably in- 
creased due to the effect of parametric resonance in " castle wall" magnetic field. Assume 
that the magnetic field is rather small. 

{^B a f = E 2 a « A 2 , (38) 

Then from the resonance condition of eq. (37) it follows that 

0! + 2 = nk, keN. (39) 

If we introduce the mean oscillation frequency as 

u = X T \ (40) 

then the resonance condition can be rewritten in the form 



n = f , (4i) 

where Q = ^ is the magnetic field variation frequency (T = Ti + T 2 ). The last formula 
again reveals the well known feature of the parametric resonance in an oscillating system: 
the resonance occurs when the doubled frequency, 2u>, is equal to the frequency of a 
parameter variation multiplied by an integer number. 

Consider the transition probability, eq.(36), in the case of the parametric resonance 
(39). Thus from (33) for this case we get 

X 2 = ±-As\El + s\E\ + 2 Sl E lS2 E 2 (-l) k ). (42) 

P " " 

Following the analysis performed in the first paper of ref. [10], one can show that the 
optimal conditions for having transition probability are realised when 

(pa = ^ + nk a , a = 1,2, (43) 

and k\ + k 2 > . Accounting (38) and (43) from (42) it follows that 



\X\ 



Ei — E 2 



< 1. (44) 



P 
Thus, for the transition probability (36) in the case of parametric resonance we get 



2 / „ ( E l ~ E l) 



P(nT) = sin 2 n v ' = sin 2 (2n(^i - 6 2 )). (45) 



In derivation of the last expression we use the conditions # a <C 1. The value of p being 

7rfc 

r 



equal to ^, the formula (45) can be rewritten in the form 



P(nT) = sin 2 {El ~ E2) t) (46) 

y i\k J 

It follows that the parametric enhancement of neutrino oscillations becomes most suffi- 
cient when k — 1. The similar effect exists for the parametric resonance in mechanical 
oscillating system that underlines the analogy mentioned above. 

Let us perform numerical estimates of possible effect of the parametric resonance 
amplification of neutrino oscillations in "castle - wall" magnetic field. We assume for 



simplicity that T\ = T2 = D and neglect matter effects. From the resonance condition 
(41) follows that 

D = |*§. (47) 

To obtain a maximal effect we chose k — 1. Then, for the particular choice of the neutrino 
characteristics Am 2 = O.leV^ 2 , E = O.leV, 9 vac = 0, we get D = \m. If the magnetic 
field profile corresponds to a chain of solenoids with opposite direction of electric currents, 
that results in B\ = —B2 = B in terms of magnetic field direction, then for the oscillation 
probability we get 

P(nT) = sin 2 (4n#), (48) 

where 29 = fiB/^-^-) For the neutrino transition magnetic moment fi = lO -10 /^ and 
the strength of magnetic field B = 10 7 G, we obtain 26 ~ 2.3 ■ 10~ 5 . Thus, the probability 
is given by 

P(nT) = sin 2 (4.6i0~ 5 n), (49) 

The probability attains the value P(nT) ^0.1 (10% of the total number of neutrinos are 
converted) when the number of periods equals to 7000. 
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